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ABSTRACT 

One of the puzzling features of solar magnetism is formation of long-living 
compact magnetic structures; such as sunspots and pores, in the highly turbu- 
lent upper layer of the solar convective zone. We use realistic radiative 3D MHD 
simulations to investigate the interaction between magnetic field and turbulent 
convection. In the simulations, a weak vertical uniform magnetic field is imposed 
in a region of fully developed granular convection; and the total magnetic flux 
through the top and bottom boundaries is kept constant. The simulation results 
reveal a process of spontaneous formation of stable magnetic structures, which 
may be a key to understanding of the magnetic self-organization on the Sun and 
formation of pores and sunspots. This process consists of two basic steps: 1) for- 
mation of small-scale filamentary magnetic structures associated with concentra- 
tions of vorticity and whirlpool-type motions, and 2) merging of these structures 
due to the vortex attraction, caused by converging downdrafts around magnetic 
concentration below the surface. In the resulting large-scale structure maintained 
by the converging plasma motions, the magnetic field strength reaches ~ 1.5 kG 
at the surface and ~ 6 kG in the interior; and the surface structure resembles 
solar pores. The magnetic structure remains stable for the whole simulation run 
of several hours with no sign of decay. 

Subject headings: sunspots — Sun: magnetic fields 

1. Introduction 

Sunspots and pores represent one of the oldest and most intriguing problem of solar 
magnetism. Despite the long history of observational and theoretical investigations (e.g. 
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Bray &: Loughhead I Il964j : iMoore fc Rabin I Il985l : iBrandenburg et al. 1 12010| ) the mechanism 
of their formation is still open. Our understanding of MHD processes on the Sun is getting 
significantly improved with the rapid progress in observational instruments, data analysis, 
methods and numerical modeling. For example, the data obtained by helioseismology have 
provided initial information about the structure and dynamics of convective flows around 



sunspots and emerging magnetic flux beneath the solar surface (e.g. iKosovichev et al.l 12000 



Zhao et al. 


2001; 


Kosovichev 


2009) 



telescopes and Hinode space mission have provided detail ed data about the fila r nentary mag- 



netic structures and flow dynamics on the surfa ce (e.g. Ilchimoto et al.l 120071 : iBonet et al. 



20081 : lAttie et al. I l2009l : iBalmaceda et al. Il2010[ ). In addition, important support for the 



understanding and interpretation of the observations is given by "realistic" radiative MHD 
numerical simulations, which are based on the first principles and take into account all essen- 
tial physical processes. The recent progress in the numerical modeling has made it possible 



to reproduce in simulat ions many observationa. 



effects in the quiet Sun region, sunspots 



and active regions 


(e.g. ; 


stein & Nord 


lund 2001: 5 


5chu 


ssler k Voder 2006: . 


Jacoutot et al. 


2008a 


b: Martfnez-Svkora et al. 2008: 


Kitiashvili et al. 


2009a, 2010), magnetic 


flux emerging 
le whole mag- 


(e.g. [ 


shibata et al. 


1989; 


Cheune et al. 


2008: Stein et al. 2009) and even tl 


netic structures, such as pores and sunspots (e.g. 


Knolker & Schiissler 198^ 


^: Stein et al. 



20031 : iBercik et al. I l2003l : iRempel et al.ll2009[ ). However, most of the modeling has been 
done by setting up the initial conditions with already existing magnetic structures, e.g. a 
horizontal flux tube for the modeling of magnetic flux emerging, or a vertical flux tube with 
strong field for the sunspot/pore structures simulations. It seems that so far only one study 
succeeded in reproducing a spontaneous formation of a micropore-like magnetic structur e 



from an initially uniform field in the turbulent convection of the Sun (jStein et al. II2003I ). 



Howeve r, the lifetime of thi s structure was rather sh ort, only "few conve ctive turnover time 



scales" ( IBercik et al. 1120031 ). Similar calculations by lVogler et al. I (120051 ) for a substantially 



shallower convective layer did not show the structure formation. 

Here, we present new results of the realistic MHD simulations that show a process of 
spontaneous formation of a stable pore-like magnetic structure from an uniform magnetic 
field, and discuss the physical mechanism of the structure formation, and its dynamics and 
evolution. 



2. Numerical setup 



For the simulations we used a 3D radiative M HD code, "SolarBox" , developed by A. 
Wray at NASA Advanced Supercomputing Division (j Jacoutot et al.|[2"o08a b: Kitiashvili et al 



- 3 - 



2009bl ). The code is built for 3D simulations of compressible fluid flows in a magnetized and 
highly stratified medium of top layers of the convective zone and the low atmosphere, in 
the rectangular geom etry. The code has been t horoughly tested and compared for test runs 



with a similar code of IStein &: Nordlund I (l200l[ ). The code is based on the Large-Eddy Sim- 



ulation (LES) approach, and solves the grid-cell averaged equations of the conservation of 
mass, momentum and energy. It takes into account the real-gas equation of state, ionization 
and excitation of all abundant species, and magnetic effects. A unique feature of the code is 
implementation of various sub-grid scale t urbulence models. For this particular simulations 



we use the minimal hyperviscosity model fjJacoutot et al.ll2008a( ) 



The simulation results are obtained for the computational domain of 6.4 x 6.4 x 5.5 Mm 
with the grid sizes: 50 x 50 x 43 km, 25 x 25 x 21.7 km and 12.5 x 12.5 x 11 km (128^ x 127, 
256^ X 253 and 512^ x 505 mesh points). The domain includes a top, 5 Mm-deep, layer of 
the convective zone and the low atmosphere. The lateral boundary conditions are periodic; 
and the top and bottom boundaries were closed for flows and maintain the constant total 
magnetic flux. The results have been verified by increasing the computational domain size to 
12.8 Mm in the horizontal directions. The initial uniform magnetic field of various strength 
1, 10 and 100 G, was superimposed over the fully developed granular convection. The 
computation runs were up to 8 hours of solar time. 

We describe first the modeling of the solar convection without magnetic field, which 
prepares the initial conditions, and then, the simulations with the superimposed weak vertical 
magnetic field. 



2.1. Simulation of the quiet Sun regions 

Figured] shows snapshots for temperature (left column) and density (right) at the surface 
for the case without magnetic field. As usual, the convective motions develop a characteristic 
granulation pattern with the relatively hot and less dense upfiowing plasma in the middle 
of the granular cells, and the lower temperature and higher density downfiowing plasma at 
the intergranulation boundaries (dark lines of granulation). An interesting feature of the 
convective flows is the formation of whirlpool-like motions of different sizes (~ 0.2 — 1 Mm) 
and lifetimes (~ 15 — 20 miij^) at the vertexes of the intergranular lanes. The vortical motions 
are particularly well seen in the density variations. The centers of the whirlpools are seen as 
dark dots (indicated by white arrows in Fig. [T)d) in the intergranular space. The evolution of 



^ The lifetime can be longer. It is difBcult to estimate it accurately because during its evolution a vortex 
can significantly change the shape and size, and sometimes almost disappear and then rise again. 
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these vortices is ultimately related to the dynamics of convective motions in the domain. The 
convective flows may sometime collect swirls in a local area, the n merge and destro y then i. 
Such vortical structures in simulations were first described by IStein &: Nordlund I ( 119981 ) . 
They showed that stronger vortices usually correlate with downflows, and this is also found 
in our results. 

From time to time, convection creates pretty big whirlpools, as the one indicated by 
square in Fig. [H which can swallow up other smaller swirls around them. The big swirls are 
usually easy to see also in the surface temperature and intensity variations. The detailed 
structure of a large whirlpool is shown in Figs [Tt-f. The whirlpool structure is characterized 
by: 1) formation of "arms" of higher density that correlates with lower temperature; 2) a 
pronounced vortical structure of the velocity flow; 3) increased magnitude of the horizontal 
velocity up to 7 - 9 km/ s; 4) a sharply decreased density in the central core of the vortex, and 
a slightly higher temperature than in the surrounding. The typical depth of large swirls is 
about 100 - 200 km. Inside the whirlpool shown in Fig. [1^, we can see a higher temperature 
tube-like structure, but it is unstable (in comparison with the whirlpool lifetime), and can 
be destroyed during the swirling motions. 

The vortic al motions in the solar g ranulation have been detected in high-resolution ob- 
servations (e.g. iPotzi fc Brandt 1 120051 ) . and the observational results generally agree with 
the simulations. In particular, recent observations of a quiet region, near the solar disk 
center detected magnetic bright points followin g a logarithmic spir al trajectory around in- 



tergranular points and engulfed by a downdraft ( iBonet et al. II2008I ). The observations were 



interpreted as vortical flows that affects the bright point motions. These whirlpools have 
the size ^0.5 Mm and the lifetime of about 5 min, wi thout preferred direction of rotation. 



The d istribution of vortices stu died from the ground (IWang et al. 



20051 . 120071 : iBonet et al. I l2008l ) and space observations ( lAttie et al. 1 120091 ) shows strong 



1995 



Potzi fc Brandt 



preferences to concentr a tion i n regions of downflows, particulary at the mesogranular scale 
( IPotzi fc Brandt I l2005l 120071 ). Our simulations for the domain of the horizontal size of 
12.8 Mm also show a tendency of concentration of vortices on a mesogranular scale. We plan 
to discuss this effect in a separate paper. Here we focus on the links between the whirlpools 
and magnetic structure formation. 



2.2. Spontaneous formation of magnetic structures 

To investigate the process of magnetic fleld structuring in the turbulent convective 
plasma we made a series of simulations for the initial vertical uniform magnetic fleld, Bzq, 
varying from 1 to 100 G, different computational grids and domain sizes. Qualitatively 
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the simulation results are very similar in all these cases, and show the formation of stable 
magnetic pore-like structures. In Figures [2] and [3] we present the results for the case of 
Bzq = 100 G, the grid size of 25 km, and the domain size of 6.4 x 6.4 x 5.5 Mm, for which 
we have made the longest run (> Shours). The periodic lateral boundary conditions allow 
us, for the illustration purpose, to shift the horizontal frame so that the structure is located 
close to the center. As we see in the simulations, the structure can be formed in any place 
of our computational domain, but usually the process starts at one of the strongest vortices. 

Figure [2] shows snapshots of the vertical magnetic field (background color image), hori- 
zontal flow field (arrows), and the vorticity magnitude (black contour lines) for four moments 
of time: 3, 10, 20 and 60 min after the moment t = 0, when the 100 G vertical field was 
uniformly distributed in the computational domain. During few minutes the magnetic field 
is swept into the intergranular lines, and is significantly amplified up to ~ 500 — 1000 G. The 
vortices and magnetic field get concentrated at some locations in the intergranular lanes, 
where they are deformed and became elongated (or elliptically shaped) along the intergran- 
ular lines. The process of formation of a large-scale magnetic structure starts at a strongest 
vortex in our domain. Since our computational domain is periodic we choose the horizontal 
coordinates in such a way that this vortex is located in the middle of Figure [2^. This large 
swirl sweeps magnetic field and also becomes stretched by strong horizontal shear flows. 
The whirlwind causes deformation of the intergranular space, and creates a cavity of low 
density, temperature and pressure. The cavity expands and increases the accumulation of 
magnetic field (Fig. [2^). A similar process of magnetic field concentration, sweeping, twist- 
ing and stretc hing by vortical mo tions in the intergranular lane was initially observed in the 



simulations of IStein et al. I (120031 ). 



During the next few minutes the deformation of the "parent" vortex continues; then 
it gets destroyed on the surface by t = 10 min (Fig. [2)d), but leaves strong downdraft 
motions in the interior (Fig. ^jp). The process of accumulation of magnetic flux in this area 
continues. The local concentrations of magnetic field and vorticity get stronger and are 
moved by convective motions in the direction of the initial cavity, into the region where the 
gas pressure remains systematically low due to the downdrafts. In Fig. |2)d several small- 
scale structures can be seen in this area. However, it is difficult to recognize the center 
of the attraction on the surface. Collisions of flows coming from different directions create 
additional vortices, and this seems to accelerate the accumulation process. As a result, the 
different small magnetic structures join together in a magnetic conglomerate that continues 
to attract other magnetic micro structures (Fig. [2t) and becomes more compact (Fig. [211). 
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This process of the magnetic structure formation is particularly well illustrated in the movi^. 
It shows a correlation between the distribution of vortices and the areas with concentrated 
magnetic field elements, as well as the mutual influence of convection, vortices and magnetic 
field. 

Figure |3] shows a vertical slice for the same moments of time. The arrows show the 
velocity field calculated from the x- and ^-components. Thus, we can see that the initially 
uniform magnetic field has local concentrations in the regions of the near-surface downdraft, 
and has a very fragmentary structure during the first few minutes (Fig. [3^). The continuing 
and extending into the deeper layers downflows are accompanied by the local concentrations 
of the magnetic field strength (Fig. Wp)i which grows after ~ 20 minutes to ~ 2 kG below 
the surface (Fig. |3t). By t = 60 min, the magnetic field is mostly concentrated in a single 
flux structure with a maximum field strength of about ~ 4 kG at the depth of ~ 1 — 
4 Mm, and ~ 1.4 — 1.5 kG at the surface (Fig. |3li). The magnetic field is weaker and more 
disperse near the bottom of our domain, which is impenetrable for flows. The fiux-tube 
interior repres ents a cluster- type structure (Figs. [2li, |3]i and H]), as initially predicted by 



Parker I ( 1l979l ) and observed on a larger scale for a sunspot by helioseismology ( IZhao et al. 
2010). In our simulations the cluster structure is represented by internal field concentrations 
(flux tubes), 100 — 200 km thick, in which the field strength reaches 6 kG after 1 hour 
(Fig. H^). The velocity distribution shows strong, often supersonic, downflows around the 
magnetic structure. Inside the magnetic structure the convective flows are suppressed by 
strong magnetic field. However, despite the weak velocities (~ 0.1 — 0.2 km/s) there are 
very small elongated convective cells resembling the umbral dots observed at the surface. 
The distribution of density fluctuations (Fig. HJd) shows the following basic properties: a) 
decrease of density inside the magnetic structure; b) fine needle-like structurization; c) a 
thin near-surface layer of slightly higher density; and d) higher density around the structure, 
particularly in the deep layers of the domain (Fig. HJd). We have followed the evolution of 
the magnetic pore-like structure for more than 8 solar hours, and did not see any indication 
of its decay. However, the shape and other properties fluctuate during this evolution. The 
decay of this structure is probably prevented by keeping the total magnetic flux constant 
during the run. We repeated the simulation when the initial uniform magnetic field was 
introduced at the different moment of time. In one case we observed the formation of two 
separate magnetic structures, which later merged together, but the whole process was very 
similar. 



See also movie: http://soi.stEinford.ed.u/~irina/fig2swirls_movie.inpeg 
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3. Conclusions 



Initi ally, the vortical mo tions in quiet- Sun regions were detected on a large scale of 
~ 5 Mm ( iBrandt et al. Ill988l ). With the developr nent of instrumenta tion it became possible 
to observe small-scale swirls in the pho tosphere (jBonet et al. 1 120081) (with size < 0.5 Mrn ) 
and also in the chromosphere ~ 1.5 Mm (jWedemeyer-Bohm &: Rouppe van der Voort II2009I ). 
Very recently, the high-resolution observations revealed a process of dragging of small-scale 
magnetic concentrations toward the center of a convective vortex motion in the photosphere 
(IBalmaceda et al. 1120101 ). 

Our simulations show that the small-scale vortices representing whirlpool-type motions 
at intersections of the intergranular lanes may play important roles in the dynamics of the 
quiet-Sun and magnetic regions. Our simulations show that the process of formation of small- 
scale magnetic structures and their accumulation into a large-scale magnetic structure is 
associated with strong vortical downdrafts developed around these structures. The resulting 
stable pore-like magnetic structure has the highest field strength of ~ 6 kG at the depth of 
1 — 4 Mm and ~ 1.5 kG at the surface. It has a cluster-like internal structurization, and 
seems to be maintained by strong downdrafts converging around this structure and extending 
into the deep layers. Our simulations show that this internal dynamics plays a critical role in 
the magnetic self-organization of solar magnetic fields and formation of large-scale magnetic 
structures. 
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Fig. 1. — Snapshots of granular convection at the surface for the simulations without mag- 
netic field, and the horizontal resolution of 12.5 km: temperature (left column) and density 
(right). The black square indicates a large whirlpool, the horizontal and vertical structure 
of which are shown in panels c) - f). Black arrows show the flow velocity. White arrows in 
panel b) point to the centers of some vortices (dark low-density points). 
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Fig. 2. — Snapshots of the surface distribution of vertical magnetic field (color background), 
horizontal flows (arrows) and vorticity magnitude (black contour lines) for four moments 
of time: 3, 10, 20 and 60 min, from the moment of initiation of a uniform magnetic field 
{Bzo = 100 G). 
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Fig. 3. — Vertical snapshots of the vertical magnetic field (color background) and the hori- 
zontal fiows (arrows) for t = 3, 10, 20 and 60 min. 
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Fig. 4. — Vertical slices though the magnetic structure at t= 84 min for the 12.5 km resolu- 
tion: a) the vertical component of magnetic field (color background), flow velocity (arrows), 
magnetic field lines (contour lines); and b) variations of density with respect to a mean den- 
sity profile of the convection simulations without magnetic field. The velocity vectors and 
magnetic field lines in panel a) are calculated from the corresponding x- and ^-components. 



